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Abstract 


Hopf  bifurcation  theory  for  an  oscillator  subject  to  a  weak  feedback  but 
a  large  delay  is  investigated  in  detail.  The  problem  is  motivated  by  semicon¬ 
ductor  laser  instabilities  which  are  generated  by  optical  feedbacks.  A  specific 
laser  problem  is  examined  for  its  Hopf  bifurcation.  Because  of  the  large  de¬ 
lay,  a  delayed  amplitude  appears  in  the  slow  time  amplitude  equation  and 
leads  to  new  bifurcations  to  periodic  and  quasiperiodic  states.  We  determine 
expressions  for  all  bifurcation  points  and  show  how  the  bifurcation  points  to 
quasiperiodic  oscillations  emerge  from  double  Hopf  bifurcation  points.  Dif¬ 
ferent  cases  of  bistability  between  steady,  periodic  and  quasiperiodic  regimes 
are  possible  which  we  study  numerically.  Finally,  the  validity  of  the  Hopf 
bifurcation  approximation  is  investigated  numerically  by  comparing  the  bi¬ 
furcation  diagrams  of  the  original  laser  equations  and  the  slow  time  amplitude 
equation. 

Keywords:  semiconductor  laser  instabilities,  system  of  delay-differential 
equations,  two-time  solution,  bifurcation  to  periodic  and  quasiperiodic  oscil¬ 
lations. 
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1  Introduction 


Semiconductor  lasers  are  the  technology  of  choice  for  many  applications  re¬ 
quiring  a  coherent  light  source  because  they  are  of  relatively  small  size,  they 
are  massively  produced  and  they  arc  easy  to  operate.  Applications  of  semi¬ 
conductor  lasers  span  a  broad  range  of  areas  from  optical  communication  to 
optical  ranging  and  sensing.  However,  semiconductor  lasers  are  extremely 
sensitive  to  optical  feedback  (OFB).  Small  amount  of  OFB  may  result  from 
the  reflection  from  an  optical  disk  or  from  the  end  of  an  optical  fiber.  This 
feedback  perturbs  the  normal  output  of  the  laser  and  generates  dynamical  in¬ 
stabilities.  These  instabilities  arc  typically  accompanied  by  unwanted  higher 
intensity  of  frequency  noise.  It  is  therefore  of  practical  interest  to  under¬ 
stand  how  the  stability  of  the  laser  depends  on  OFB.  A  simple  model  has 
been  formulated  in  1980  by  Lang  and  Kobayashi  (LK)  and  consists  of  two 
nonlinear  ordinary  delay-differential  equations  for  the  complex  electrical  field 
Y(t)  and  the  carrier  number  Z(t)  [?].  Time  t  is  measured  in  units  of  the  pho¬ 
ton  lifetime  rp.  A  term  proportional  to  Y(t  —  r)  appears  in  the  field  equation 
and  describes  the  effect  of  OFB  after  the  reffection  of  light  from  a  mirror. 
t  is  proportional  to  the  round-trip  time  from  laser  to  mirror  and  back  to 
laser.  Computer  simulations  have  shown  that  the  LK  delay-differential  equa¬ 
tions  correctly  describe  the  dominant  effects  observed  experimentally.  This 
includes  the  occurrence  of  mode  hopping  [?,  ?],  coexisting  time-periodic  at¬ 
tractors  [?,  ?,  ?]  and  different  forms  of  chaotic  attractors  called  low-frequency 
fluctuations  [?,  ?,  ?].  More  recently,  a  numerical  bifurcation  study  suggests 
that  LK  instabilities  may  result  from  cascading  bifurcations  starting  from  a 
Hopf  bifurcation  [?].  Ideally,  LK  equations  models  an  experimental  set-up 
consisting  of  a  laser  exposed  to  weak  reflection  of  light  from  a  mirror  located 
at  a  sufficiently  large  distance  to  avoid  multiple  reflections  (one  to  two  me¬ 
ters).  This  implies  a  round-trip  time  r'  from  laser  to  mirror  and  back  to  laser 
of  the  order  of  1  ns  which  is  relatively  large  compared  to  the  photon  lifetime 
tp  ~  1  ps  (i.e.,  a  dimensionless  delay  r  =  t' /tp  —  O(103)).  It  is  this  large 
delay  r  which  is  causing  the  laser  instabilities  even  for  low  feedback  rates  and 
which  motivates  our  mathematical  interest  for  the  Hopf  bifurcation  subject 
to  a  large  delay. 

Periodic  solutions  of  delay-differential  equations  and  their  behavior  for 
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large  delay  have  been  studied  mathematically  for  equations  of  the  form 


uf  —  f{u,u(t-  t),A)  (1) 

where  u  €  R.  Applications  appear  in  the  biological  sciences  [?],  chemical 
problems  [?]  and  in  nonlinear  optics  [?,  ?].  Most  of  the  solutions  computed 
numerically  correspond  to  square-wave  oscillations  exhibiting  a  period  close 
to  r  as  r  — >  oo.  The  condition  for  the  bifurcation  to  these  square-wave 
periodic  regimes  has  been  examined  mathematically  in  [?,  ?].  By  contrast 
to  these  square-wave  oscillations,  the  oscillatory  solution  of  LK  equations 
for  low  feedback  is  nearly  harmonic  in  time  [?,  ?]  and  the  frequency  of  the 
oscillations  is  0(1)  as  r  oo.  The  effect  of  a  large  r  appears  through 
a  secondary  bifurcation  to  quasipcriodic  oscillations  which  exhibits  a  low 
frequency  proportional  to  r_1.  This  typical  two-time  behavior  of  the  solution 
of  LK  problem  for  low  feedback  rates  but  large  delays  motivates  the  study 
of  specific  Hopf  bifurcation  problems  of  the  form 

v!  =  f(u,e2u(t  -  t),X)  (2) 

where  u  E  Rn  (n  >  2).  In  Eq.  (??),  e  and  r  represent  a  small  parameter 
and  a  large  parameter,  respectively.  General  Hopf  bifurcation  theories  for 
systems  of  delay-differential  equations  arc  difficult  (see  [?])  and  we  shall 
take  advantage  of  the  small  parameter  e.  Specifically,  we  assume  that  the 
e  —  0  problem  admits  a  Hopf  bifurcation  at  A  —  A0  from  u  —  0  to  a  stable 
time-periodic  state.  Our  objective  is  to  investigate  this  bifurcation  for  small 
e  and  progressively  larger  values  of  r.  To  this  end,  we  may  apply  a  two- 
time  perturbation  method  introducing  a  fast  time  t  and  a  slow  time  u=  e2t 
[?,  ?].  The  two  independent  time  variables  imply  that  u(t)  —  u(t,u )  and 
u(t  —  t)  —  u(t  —  r,  v  —  e2r).  If  r  is  0(1),  u(t  —  t,u  —  e2r)  simplifies  as 

u(t  —  r,  v  —  e2r)  =  u(t  —  r,  u)  +  0(e2)  (3) 

and  the  effect  of  the  delay  on  the  slow  time  u  will  not  appear  in  the  leading 
order  amplitude  equation.  On  the  other  hand,  if  r  is  0(e~2)  large,  we  cannot 
expand  u(t  —  t,u  —  e2r)  for  small  e  and  the  slow  time  delay  given  by  e2r  will 
appear  explicitly  in  the  bifurcation  equation  through  a  delayed  amplitude. 
As  we  shall  demonstrate,  this  slow  time  delay  is  responsible  for  multiple 
branching  of  solutions. 
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In  this  paper,  we  concentrate  on  a  simple  laser  problem  which  is  modeled 
by  equations  of  the  form  of  (??).  This  will  allow  us  to  discuss  the  validity 
of  Hopf  local  theory  by  comparing  the  numerical  bifurcation  diagrams  ob¬ 
tained  from  the  original  equations  and  from  the  slow-time  amplitude  equa¬ 
tion.  Specifically,  we  shall  consider  a  semiconductor  laser  subject  to  a  delayed 
feedback  controlled  opto-electronically  [?].  The  problem  is  mathematically 
simpler  than  LK  problem  because  the  phase  of  the  laser  held  docs  not  play  a 
role  and  the  dynamical  variables  reduce  to  the  intensity  of  the  laser  held  and 
the  electronic  carrier  density.  The  laser  equations  arc  formulated  in  dimen¬ 
sionless  form  in  Section  2  and  the  slow- time  amplitude  equation  is  derived 
in  Section  3.  In  Section  4,  we  determine  its  steady  state  solutions  and  for¬ 
mulate  stability  conditions  which  we  investigate  numerically.  We  conclude 
that  multiple  branches  of  periodic  and  quasipcriodic  states  arc  possible  and 
are  connected  at  double  Hopf  bifurcation  points.  In  Section  5,  we  evaluate 
the  validity  of  our  Hopf  bifurcation  analysis  by  comparing  numerical  bifurca¬ 
tion  diagrams.  We  also  identify  different  forms  of  bistability  between  steady, 
periodic  and  quasiperiodic  regimes  which  are  of  particular  interest  for  exper¬ 
iments.  Finally,  Section  6  discusses  the  physical  relevance  of  our  results  for 
semiconductor  lasers  instabilities. 


2  Formulation 

The  semiconductor  laser  system  and  its  optoelectronic  feedback  is  sketched  in 
Figure  1.  Part  of  the  laser  output  is  detected  with  a  high-speed  photodetector. 
The  detector  photocurrent  component  is  selected  with  a  T  bias,  amplified 
and  added  to  the  DC  pump  current  [?].  In  dimensionless  form,  the  laser  rate 
equations  modeling  this  system  are  given  by  [?] 

dY 

—  =  (1  +  ia)ZY,  (4) 

T^^P  +  i,  |  Y(t  -r)f-Z-(l+  2  Z)  \Yf  (5) 

where  Y  and  Z  represent  the  complex  electrical  held  and  the  electronic  carrier 
density.  Time  t  is  measured  in  units  of  the  photon  lifetime  rp  ( t  =  t'/rp).  a 
is  defined  as  the  linewidth  enhancement  factor,  T  =  rn/rp  is  a  ratio  of  two 
time  constants  where  rn  is  defined  as  the  carrier  lifetime,  r  =  r' /rp  is  the 
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dimensionless  delay  of  the  feedback  where  r'  is  the  cavity  transit  time  and  P 
is  the  pump  parameter  above  threshold  .  These  equations  arc  the  traditional 
semiconductor  laser  rate  equations  with  P  now  replaced  by  P  +  77  \Y(t  —  r)  |2 
which  models  the  effect  of  the  DC  coupled  optoelectronic  feedback.  Typical 
values  of  the  fixed  parameters  arc  a  ~  3  —  6,  rp  ~  2  ps,  rn  ~  2  ns  and  t 
~  1  —  10  ns  which  mean  large,  O(103)  values,  for  both  T  and  t. 

We  shall  take  advantage  of  these  two  large  parameters  in  order  to  deter¬ 
mine  an  approximation  of  the  solution  of  Eq.  (??)  and  (??).  Introducing  the 
decomposition  Y  —  \fl  exp (iip)  allow  us  to  reduce  Eqs.  (??)  and  (??)  to  two 
coupled  equations  for  the  intensity  /  and  the  carrier  density  Z.  ip  is  passively 
related  to  I  and  Z.  Furthermore,  the  large  T  parameter  which  multiplies 
the  left  hand  side  of  the  carrier  equation  (??)  can  be  removed  by  a  change 
of  variables.  Introducing  the  new  variables  x ,  y  and  s  defined  by 

I  =  P{  1  +  y),  Z  =  ^x,  s  =  ut,  (6) 


where 

w  =  fip/T  (7) 

is  known  as  the  laser  relaxation  oscillation  frequency,  the  two  equations  for 
I  and  Z  become 


/  2  P  \ 

rf  =  -V~  &  [1  +  1  +  2pyj  +rj[  1  +  y{s  -  0)] ,  (8) 

y'  =  { 1  +  y)%  (9) 

where  prime  means  differentiation  with  respect  to  s.  The  new  variables  x  and 
y  represent  the  deviations  of  I  and  Z  from  their  steady  state  values  I  —  P 
and  Z  —  0,  respectively.  The  new  parameters  £  and  9  are  defined  by 

f  XpP^  and  9  ~  UT' 

^  is  known  as  the  damping  rate  of  the  2it/lu  laser  intensity  oscillations  and  9 
is  a  rescaled  delay.  The  large  T  parameter  now  appears  in  the  small  param¬ 
eter  £  (£  is  proportional  to  u>  which  is  0{T~1/2)  small).  Eqs.  (??)  and  (??) 
are  in  an  appropriate  form  for  an  asymptotic  analysis.  These  equations  arc 
equivalent  to  the  equations  studied  in  [?,  ?]  for  a  variety  of  different  opto¬ 
electronic  devices.  In  [?,  ?],  the  objective  was  to  determine  strongly  pulsating 
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periodic  solutions  for  moderate  values  of  the  delay.  Here,  we  determine  peri¬ 
odic  solutions  which  are  nearly  harmonic  for  large  values  of  the  delay.  Recent 
attempts  to  investigate  the  large  delay  limit  for  similar  laser  equations  arc 
given  in  [?,  ?].  Simplified  amplitude  equations  were  derived  but  their  numer¬ 
ical  validity  as  providing  approximations  of  the  solutions  of  the  original  laser 
equations  have  not  been  examined.  In  this  paper,  we  determine  analytical 
approximations  of  several  branches  of  periodic  solutions  and  their  validity  is 
investigated  in  detail. 


3  Hopf  bifurcation 

In  this  section,  we  determine  a  slow  time  amplitude  equation  by  a  two-time 
analysis  of  Eqs.  (??)  and  (??)  based  on  the  limit  £  and  rj  small  but  9  large. 
The  perturbation  analysis  is  quite  simple  and  we  summarize  the  main  results. 


We  first  introduce  a  small  parameter  e  defined  by 

(  =  (ii) 

and  assume  the  following  scalings  for  the  parameters  r/  and  9 

77(e)  -  e2C,  9  =  W20.  (12) 

We  seek  a  two-time  solution  of  the  form 

x(s,  u,  e)  —  exi (s,  u)  +  e2x2{s,  v)  +  ...  (13) 

y(s,  v,  e)  =  eyi(s,  v)  +  e2y2(s,  u)  +  ...  (14) 

where  r  is  a  slow  time  variable  defined  by 

//  =  e2s.  (15) 


Inserting  (??)-(??)  into  Eqs.  (??)  and  (??),  using  the  chain  rule  d/ds  — 
d/ds  +  €ldjdv  and  noting  that  y{s  —  9)  —  y(s  —  9,  u  —  0),  lead  to  a  sequence 
of  problems  for  the  unknown  coefficients  in  (??)  and  (??).  The  solution  of 
each  problem  is  easily  determined.  We  find  that 

x  —  t  (AM  expRs)  +  c.c .)  +  0(e2),  ,  , 

y  =  —t  (iA(v)  exp(*s)  +  c.c.)  +  0(e2)  (16) 
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where  c.c.  means  complex  conjugate  and  the  slow-time  amplitude  A  satisfies 
the  following  ordinary  differential  equation 


A'  = 


-A2  A* 


A  —  iCA(v  —  0)  exp(— i9) 


(17) 


The  main  difference  with  the  case  9  —  0(1)  is  the  presence  of  the  delayed 
amplitude  A  (is  —  0)  in  Eq.  (??)  instead  of  A(is).  It  is  instructive  to  review 
the  significance  of  each  term  in  the  right  hand  side  of  Eq.  (??).  The  first  two 
terms  in  the  right  hand  side  of  Eq.  (??)  represent  the  linear  and  nonlinear 
correction  to  the  frequency  and  could  be  anticipated  from  a  study  of  the  laser 
equations  (??)  and  (??)  without  damping  (  £  =  0)  and  without  feedback 
(rj  —  0).  The  third  term  in  the  right  hand  side  of  Eq.  (??)  comes  from  the 
coefficient  multiplying  £  in  Eq.  (??).  It  represents  the  physical  damping  of 
the  laser  oscillations.  Finally  the  last  term  in  the  right  hand  side  of  Eq.  (??) 
is  the  contribution  of  the  feedback. 

Eq.  (??)  is  the  starting  point  of  our  bifurcation  analysis.  For  clarity,  the 
different  parameters  introduced  in  the  last  two  sections  are  shown  in  Table 
1.  _ 


Symbol 

Definition 

Value 

u> 

7? 

0.03 

relax,  oscill.  frequency 

e2  =  f 

,  ,  1+21*  _  L+2P 

W  IP  s/2PT 

0.06 

damping  rate 

9 

LOT  —  \j‘^rT 

31.62  to  63.25 

scaled  delay 

0 

(N  r 

4-^ 

■ — i 

II 

2  to  4 

scaled  delay 

C 

r,/e2 

—  1.6  to  1.6 

scaled  feedback  rate 

Table  1.  Laser  paramerers  and  their  typical  values.  T  —  1000,  P  —  0.5  and 
t  =  1000  to2000. 


4  Analysis  of  the  slow  time  amplitude  equa¬ 
tion 


In  this  section,  we  determine  the  steady  state  solutions  of  Eq.  (??)  and 
investigate  their  linear  stability  properties.  The  validity  of  Eq.  (??)  as  an 
approximation  of  the  solutions  of  the  original  laser  equations  (??)  and  (??) 
will  be  investigated  numerically  in  Section  5.  The  main  advantage  of  Eq. 
(??)  is  the  fact  that  periodic  and  quasiperiodic  solutions  of  the  original  laser 
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equations  now  correspond  to  steady  and  periodic  solutions  for  the  amplitude 
|A|.  It  is  mathematically  convenient  to  introduce  the  decomposition  A  — 
Rexp(i(f>)  into  Eq.  (??).  We  then  obtain  the  following  equations  for  R  and  0 


R!  =  -  [-R  +  CR(is  -  0)  sin  (-9  +  cf>(is  -  0)  -  0)] 

Li 


^'=2 


CR  -  CR(u  -  0)  cos  (-9  +  <t>(v  -  0)  -  0)  -  l-R 3 

3 


(18) 

(19) 


Before  we  analyze  these  equations  in  detail,  it  is  interesting  to  determine  the 
solution  of  the  regular  problem  with  9  —  0(1).  Equivalently,  we  solve  Eqs. 
(??)  and  (??)  with  0  —  0.  A  periodic  solution  of  the  original  laser  equations 
corresponds  to  R  —  cst  and  0  —  (f)(0)  +  Bis.  From  Eqs.  (??)  and  (??)  with 


0  =  0,  we  find  that 


C  =  C°H  = 


1 

sin(0) 


(20) 


and 

R  =  ^6(fl„  -  B) 


(21) 


where  Z?0  is  defined  by 

Ba  =  —  ^tan(^)  (22) 

and  B  <  B()  is  arbitrary.  C  —  C\  is  a  Hopf’  bifurcation  point  and  the  bifurca¬ 
tion  is  vertical  at  this  order  of  the  analysis.  A  higher  order  analysis  will  lead 
to  the  direction  of  bifurcation  but  we  do  not  need  this  result  (see  Appendix 
in  [?]).  Our  main  observation  is  that  the  Hopf  bifurcation  is  unique  if  the 
delay  9  is  an  0(1)  quantity.  As  we  shall  now  demonstrate,  increasing  9  will 
bend  the  vertical  bifurcation  branch  (??)-(??)  and  produce  new  (bifurcation 
and  isolated)  branches  of  periodic  solutions. 


4.1  Periodic  solutions 


We  consider  Eqs.  (??)  and  (??)  with  0^0  and  seek  a  solution  of  the  form 
R  —  cst  and  (f>  —  0(0)  +  Bis.  This  solution  corresponds  to  a  periodic  state 
of  the  original  laser  equations  which  exhibits  a  frequency  equals  to  1  +  e2B. 
The  exact  solution  for  R2  —  R2(C)  and  B  —  B(C)  can  be  formulated  in 


parametric  form  as 


C  — 


1 

sin  (9  +  BB)  ’ 


(23) 
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R2  —  3  [-2B  —  tan  ((9  +  BQ)/2)\  >  0  (24) 

where  the  frequency  D  is  the  parameter.  Changing  D  from  —  oc  to  +oo,  we 
obtain  C  and  R  using  (??)  and  (??).  See  Figure  2.  The  figure  show  that 
there  exist  several  branches  of  solutions.  The  limit  as  0  — »  0  is  discussed  in 
Appendix  A. 

From  (??)  and  (??)  with  R  —  0,  we  find  that  the  Hopf  bifurcation  points 
( B,C )  —  ( Bh,Ch )  satisfy  Eq.  (??)  and 

2D  +  tan  ((9  +  BQ)/ 2)  =  0.  (25) 


Thus,  we  determine  D  —  BH  by  solving  Eq.  (??)  and  then  obtain  C  —  Ch 
using  (??).  A  different  and  more  physically  interesting  way  to  investigate 
the  roots  of  Eq.  (??)  as  a  function  of  9  is  to  study  the  implicit  expression 


given  by 


9  - 


2 

1  +B$ 


TO7T 


arctan(2Z?)] 


(26) 


where  m  —  0, 1,  2...  Having  9  —  9(B),  we  determine  C  —  Ch  by  (??).  Thus, 
we  obtain  9  and  Ch  by  changing  parameter  B  continuously  from  — oo  to 
oo.  See  Figure  3.  For  each  m  ,  we  have  a  line  of  bifurcation  points  which 
corresponds  to  a  different  branch  of  periodic  solutions. 

The  direction  of  bifurcation  can  be  determined  from  an  analysis  of  Eqs. 
(??)  and  (??)  for  ( B.  C )  close  to  ( BH,CH )•  We  find  that  the  bifurcation  is 
defined  for  C  >  Ch  if 


cos (9  -\-  Bh&)  <  0,  (27) 

or  equivalently,  for  C  <  Ch  if  cos (9  +  BHQ)  >  0.  The  direction  of  bifurcation 
changes  if  cos(#  +  BhO )  —  0  which  implies,  using  (??),  that  Ch  —  1-  Since 
Ch  —  1  at  the  minimum  of  the  curve  Ch  —  Ch  (9),  the  change  of  direction 
occurs  exactly  at  that  point.  The  successive  minima  are  located  at  9  —  9rn. 
If  Ch  >  0,  we  have  verified  that  the  bifurcation  is  supercritical  (C  >  Ch)  if 
9  <  9rn  and  is  subcritical  ( C  <CH)  if  9  >  9rn. 

Figure  3a  also  shows  that  distinct  Hopf  bifurcation  lines  may  intersect. 
These  points  correspond  to  degenerate  Hopf  bifurcation  points  which  ex¬ 
hibit  two  distinct  periodic  modes  (double  Hopf  bifurcation  points).  From 
these  points,  the  solution  is  in  first  approximation  a  linear  combination  of 
the  two  periodic  modes  and  multiple  branching  of  solutions  arc  possible  at 
these  points.  We  anticipate  that  in  addition  to  the  two  pure  mode  bifurca¬ 
tions,  mixed  mode  bifurcations  are  possible.  We  verify  this  numerically  and 
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analytically  in  the  next  subsection.  This  mixed  mode  regime  is  typically 
quasiperiodic  exhibiting  the  two  pure  mode  frequencies. 

By  contrast  to  the  Hopf  points,  the  limit  points  ( B,C )  —  ( Bl,Cl ) 
lead  to  isolated  branches  of  periodic  solutions.  They  satisfy  the  condition 
sin  (9  +  BlO )  —  ±1  and  they  are  located  at 


Cl  —  ±1  and  BL 


1 

0 


(=F(1  +  4u)/t/2  -  9) 


(28) 


where  n  —  0, 1, 2, ... 

4.2  Quasiperiodic  bifurcations 

We  next  investigate  the  stability  of  the  periodic  states  given  by  (??)  and 
(??).  A  detailed  stability  analysis  of  these  solutions  is  difficult  because  the 
characteristic  equation  is  transcendental.  But  we  may  examine  the  case  of 
small  eigenvalues  and  the  case  of  purely  imaginary  eigenvalues. 

From  the  linearized  equations,  we  determine  the  characteristic  equation 
for  the  growth  rate  a.  Assuming  a  small,  we  then  obtain  an  expression  for  a 
of  the  form 

a  ~  acos(9  +  BO)R 2  (29) 

where  a  is  always  positive,  if  all  the  remaining  eigenvalues  have  a  negative 
real  part,  (??)  implies  that  the  Hopf  bifurcation  branch  is  stable  near  its 
bifurcation  point  C  —  Ch  (limit  R  — >  0,  cos [9  +  BO)  fixed)  provided  that 
(??)  is  satisfied,  i.e.,  for  a  supercritical  Hopf  bifurcation.  The  expression  (??) 
also  implies  the  local  stability  of  one  of  the  two  branches  of  periodic  solutions 
that  emerges  from  a  limit  point  C  —  CL  (recall  that  cos(0  +  BO)  —  0  at 
C  —  Cl  and  consider  the  limit  |cos(#  +  BO) \  — >  0,  R  fixed  of  Eq.  (??)). 

We  next  examine  the  conditions  for  a  Hopf  bifurcation  from  the  R  —  c-st 
solution  meaning  a  secondary  bifurcation  to  quasiperiodic  oscillations  for 
the  laser  intensity.  Substituting  a  —  into  the  characteristic  equation,  we 
obtain  two  conditions  for  the  critical  feedback  rate  C  and  frequency  //,.  They 
are  given  by 

— 4^i2  +  C 2  ((cos {jiO)  —  l)2  —  sin2(p,0)) 

+(cos(Ai0)  -  1)C cos(#  +  BO)\Rz  (30) 

—  sin(//,0)4/i  —  0 


and 
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(31) 


— 2C'2(cos(p0)  —  1)  sin (//,0) 

—4 /j,  (cos(p0)  —  1) 

—  sin(/i0)(7  cos  (9  +  i?0)|i?2  —  0 

where  R  and  D  are  related  to  C  by  (??)  and  (??).  The  solution  of  these 
equations  has  been  determined  numerically.  See  Figure  4.  We  find  successive 
lines  of  quasiperiodic  bifurcation  points  that  all  emerge  from  the  double  Hopf 
bifurcation  points.  These  points  arc  denoted  by  P\  and  P4  in  Figure  4.  I11 
Appendix  B,  we  show  analytically  that  these  lines  emerge  from  the  double 
Hopf  bifurcation  points  with  the  new  frequency  fi  —  B2  —  B]_.  and  B2 
denote  the  two  Hopf  frequencies  at  the  double  Hopf  point.  This  is  consistent 
with  a  bifurcation  theory  near  a  multiple  eigenvalue  since  we  expect  a  solution 
of  the  form 


A~  Ri  exp  (i(Biis  +  0X))  +  R2  exp  (i(B2v  +  02))  (32) 

near  this  point.  At  the  end  of  Appendix  B,  we  show  that  the  secondary 
bifurcating  solution  has  indeed  the  form  of  (??).  A  detailed  analysis  of  the 
solutions  near  the  double  Hopf  bifurcation  point  is  beyond  the  scope  of  this 
paper  devoted  to  a  single  Hopf  bifurcation  subject  to  a  large  delay.  I11  the  next 
section,  we  investigate  numerically  the  bifurcation  possibilities  suggested  by 
Figure  4  and  emphasizes  cases  of  coexisting  stable  states  (bistability). 

5  Numerical  Bifurcation  diagrams. 

The  objectives  of  this  section  are  twofold.  First,  we  investigate  numerically 
several  cases  of  coexisting  (steady  or  not)  attractors.  These  possibilities  are 
of  physical  interest  because  they  can  be  observed  experimentally.  Second,  we 
examine  the  validity  of  our  approximations  by  comparing  the  numerical  bi¬ 
furcation  diagram  of  the  original  laser  equations  and  the  bifurcation  diagram 
obtained  from  the  amplitude  equation. 

We  first  determine  the  bifurcation  diagram  of  the  solutions  of  Eq.  (??)  for 
different  values  of  9.  The  main  advantage  of  this  equation  is  that  steady  and 
periodic  states  correspond  to  periodic  and  quasiperiodic  states  of  the  original 
laser  equations.  The  simplicity  of  Eq.  (??)  allowed  us  to  have  analytical  ex¬ 
pressions  for  all  the  steady  state  branches.  We  shall  limit  our  analysis  to  the 
range  60  <  9  <  70  and  consider  C  >  0.  We  follow  the  predictions  of  Figure 


12 


4  and  our  results  are  summarized  by  the  three  qualitative  diagrams  shown 
in  Figure  5  and  by  the  detailed  diagram  given  in  Figure  6.  In  Figure  5,  we 
concentrate  on  the  gradual  change  of  the  first  Hopf  bifurcation  branch  as  the 
delay  9  is  progressively  increased.  Figure  6  examines  a  case  of  two  coexisting 
branches  of  periodic  solutions.  All  our  diagrams  arc  based  on  the  numerical 
determination  of  the  stable  solutions  of  the  slow  time  amplitude  equations 
(??)  and  (??)  as  well  as  the  exact  analytical  expressions  of  the  branches  of 
periodic  solutions  given  by  (??)  and  (??).  In  Figure  6,  we  determine  the 
stable  solutions  from  the  amplitude  equations  as  well  as  from  the  original 
laser  equations.  Four  key  points  in  Figure  4  mark  qualitative  changes  of  the 
bifurcation  diagrams.  As  described  in  Section  4.2,  P\  and  P±  are  double 
Hopf  bifurcation  points  where  lines  of  secondary  bifurcation  points  appear. 
P-2  refers  to  the  change  of  direction  of  the  Hopf  bifurcation  and  was  ana¬ 
lyzed  in  Section  4.1.  it  verifies  the  condition  Cp  —  Cqp  but  corresponds  to 
distinct  bifurcation  points  (different  amplitudes  of  the  solution). 

The  simplest  case  occurs  if  0Fl  <  9  <  9p2  (Figure  5a).  A  supercritical 
Hopf  bifurcation  is  followed  by  a  secondary  bifurcation  to  quasiperiodic  os¬ 
cillations.  If  9p2  <  0  <  9p:s  (Figure  5b),  the  sequence  of  bifurcations  is 
similar  to  the  sequence  shown  in  Figure  5a  except  that  the  Hopf  bifurca¬ 
tion  is  subcritical  and  allows  the  coexistence  of  a  stable  steady  state  and 
a  stable  periodic  regime.  This  coexistence  becomes  richer  if  9ps  <  9  <  0pA 
(Figure  5c):  in  addition  to  coexisting  steady  and  periodic  states,  we  may 
now  have  coexisting  stable  steady  and  quasiperiodic  regimes.  If  9  >  9 piy 
new  branches  of  periodic  and  quasiperiodic  solutions  appear.  An  example  is 
studied  quantitatively  in  Figure  6. 

Specifically,  we  compare  the  bifurcation  diagrams  of  the  periodic  and 
quasiperiodic  states  obtained  from  the  full  laser  equations  (??)  and  (??) 
(Figure  6a)  and  from  the  slow  time  amplitude  equations  (??)  and  (??)  (Figure 
6b).  In  Figure  6b,  we  observe  a  stable  Hopf  bifurcation  branch  that  emerges 
from  zero  at  r/  ~  0.076.  This  bifurcation  is  followed  by  a  secondary  bifurcation 
to  a  branch  of  stable  quasiperiodic  oscillations  at  rj  ~  0.086.  Near  rj  ~  0.1, 
we  note  a  tertiary  bifurcation  to  a  more  complicated  state.  In  addition  to  the 
Hopf  bifurcation  branch,  there  exists  an  isolated  branch  of  periodic  solutions 
that  appears  from  a  limit  point  located  near  rj  ~  0.064.  This  branch  then 
exhibits  a  secondary  bifurcation  to  a  quasiperiodic  state  at  rj  ~  0.068  which 
is  itself  followed  by  a  tertiary  bifurcation  to  more  complex  quasiperiodic 
oscillations  ( rj  ~  0.082).  We  compare  this  bifurcation  diagram  to  the  diagram 
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obtained  by  solving  directly  the  laser  equations  (Figure  6a).  The  agreement 
is  excellent  for  the  low  amplitude  branches.  The  agreement  becomes  more 
qualitative  for  the  higher  amplitude  branches  of  periodic  solutions,  as  we 
may  expect,  but  the  same  sequence  of  bifurcation  transitions  is  observed 
in  both  diagrams.  Note  that  e  =  0.17  for  our  parameter  values.  We  have 
verified  that  the  numerical  agreement  between  diagrams  improves  for  lower 
e  (higher  T). 


6  Discussion 


The  main  objective  of  this  paper  was  to  evaluate  the  validity  of  Hopf  bifur¬ 
cation  theory  for  a  system  of  two  delay-differential  equations  modeling  semi¬ 
conductor  lasers  subject  to  a  weak  optical  feedback  but  exhibiting  a  large 
delay.  To  this  end,  we  have  compared  the  numerical  bifurcation  diagrams 
obtained  from  the  original  laser  equations  and  from  the  slow  time  amplitude 
equations.  In  addition,  we  have  determined  useful  analytical  expressions 
for  branches  of  periodic  solutions  and  formulate  conditions  for  their  bifurca¬ 
tions.  These  expressions  considerably  helped  the  numerical  study.  Previous 
attempts  [?,  ?]  have  tried  to  take  advantage  of  this  large  delay  in  order  to 
derive  simplified  laser  equations  but  have  never  estimated  the  asymptotic 
validity  of  their  approximations. 

An  important  result  of  our  bifurcation  analysis  is  the  observation  that 
a  large  delay  may  generate  multiple  branches  of  periodic  and  quasiperiodic 
oscillations  even  for  low  feedback.  These  multiple  Hopf  or  isolated  branching 
of  solutions  are  predicted  by  our  amplitude  equation.  Because  the  derivation 
of  this  amplitude  equation  is  relatively  easy,  Hopf  local  theory  was  an  useful 
guide  for  our  numerical  study  of  the  laser  equations. 

The  multiple  branching  of  solutions  docs  not  necessarily  lead  to  a  chaotic 
laser.  We  have  shown  that  different  forms  of  bistable  responses  are  possible 
near  the  first  instability  of  the  laser.  These  cases  are  of  practical  interest 
for  semiconductor  laser  experiments  which  arc  mainly  based  on  the  deter¬ 
mination  of  power  spectra.  Jumps  between  steady,  periodic  or  quasiperiodic 
regimes  can  be  observed  by  the  changes  of  frequencies. 

We  have  shown  that  a  secondary  bifurcation  to  stable  quasiperiodic  os¬ 
cillations  results  from  the  interaction  of  two  successive  primary  Hopf  bifur¬ 
cations.  The  new  frequency  equals  the  difference  between  the  two  individual 
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Hopf  frequencies.  Our  analysis  encourages  investigating  the  possible  sec¬ 
ondary  bifurcations  in  LK  equations.  For  LK  problem,  it  is  suspected  that  a 
secondary  bifurcation  to  unstable  quasiperiodic  oscillations  may  explain  the 
coexistence  of  two  stable  primary  Hopf  bifurcation  branches  [?]  [?]. 
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8  Appendix  A.  The  small  delay  limit  (0—^0). 

In  this  appendix,  we  examine  the  behavior  of  the  Hopf  bifurcation  and  limit 
points  for  small  0  using  Eqs.  (??)  and  (??).  As  0— >  0,  all  Hopf  bifurca¬ 
tion  branches  moves  to  higher  values  of  D  (B  —  Q(0-1))  except  one  which 
converges  to  the  vertical  branch  (??).  More  precisely,  we  find  that 


where  the  new  bifurcation  point  Ch  —  Cjj  is  shifted  as 

CW  (0)  =  C°H{1-  BB0  cot(9)) .  (34) 

Thus  the  bifurcation  becomes  vertical  as  0  0  (C  —  C #  —  0(0)  if  R  — 

0(1)).  All  the  other  Hopf  bifurcation  branches  moves  to  infinity  as  ©  — >  0. 
We  find 

CH  -  ^  ((1  +  2m) 7r  -  9)  (35) 

R  =  x/6  (c  ~  CH(m)).  (36) 


9  Appendix  B.  Double  Hopf  bifurcation  points. 

For  critical  values  of  the  delay  9  (points  Pi  and  in  Figure  4),  a  Hopf 
bifurcation  point  may  exhibits  two  periodic  modes  with  the  distinct  frequen¬ 
cies  B  —  Bi  and  B  —  B2.  From  Eq.  (??),  we  find  the  simple  condition 
sin(0  +  .010)  —  sin(0  +  B20)  which  implies  that  either 

(B2  -  Bl)(-)  -  2rrm  (37) 

where  rn  —  1,2,..  or 


29  +  (B2  +  B  i)0  =  7r  +  4n7r 


(38) 


where  n  —  0,1,2,...  We  first  examine  Case  (??).  From  Eq.  (??),  we  know 
that  Bi  and  B2  verify  the  equations 


2B\  +  tan  ((9  +  Z?40) /2)  —  0 


(39) 
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(40) 


2B2  +  tan  {{9  I  -  B2Q)/2)  —  0. 
Substracting  (??)  and  (??),  we  obtain 


2 {By  -  B2)  + 


sin  {{By  -  B2)&/ 2) 
cos  (( 9  +  B !0) /2)  cos  ((0  +  B20) /2) 


Using  now  (??),  we  find  that  Eq.  (??)  cannot  be  satisfied.  We  next  consider 
Case  (??).  We  may  rewrite  the  denominator  in  (??)  as  a  sum  of  two  cosine 
functions  and,  using  (??),  we  obtain 

{By  -  B2)  +  tan  {{By  -  B2)S/ 2)  =  0.  (42) 

Eq.  (??)  is  an  equation  for  the  difference  between  the  two  frequencies.  Sim¬ 
ilarly,  we  obtain  an  equation  for  the  sum  of  the  two  frequencies  by  adding 
(??)  and  (??).  Using  (??)  again,  we  find 


7r  +  Arm  —  29 

0 


COS  {{By  I?2)0/2) 


Introducing  z  =  {By  —  B2)0/ 2  and  eliminating  9  and  0  —  (30  using  (??)  we 
obtain  an  equation  for  z  only: 

P  [(7T  +  4n7r)  sin(^)  —  2 z\  +  Az  cos(^)  —  0.  (44) 

Having  z  from  Eq.  (??),  we  determine  $  or  0  —  $/£  by  using  (??): 


^  -  ^) ' 

For  example,  with  P  —  0.5  and  T  —  1000,  we  find  z  ~  2.28,  9  ~  61.4 
(equivalently,  0  ~  3.88)  which  gives  By  —  B2  ~  1.17.  This  point  corresponds 
to  point  Py  in  Figure  4. 

We  next  consider  conditions  (??)  and  (??)  for  the  secondary  bifurcation 
point  and  show  that  the  frequency  of  the  quasiperiodic  oscillations  is  p  = 
{By  —  B2)  at  the  double  Hopf  bifurcation  point.  To  this  end,  we  consider 
Eqs.  (??)  and  (??)  evaluated  at  this  point  {R  —  0  and  C  —  C'h)  ■ 

f-P  +  Cjj  sin2  (p.0/2)  cos (p0)  +  sin(//,0)/r  =  0  (46) 
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sin2(p0/2)  ( ' C\  sin(//,0)  +  2 fij  =0  (47) 

Eq.  (??)  is  satisfied  if  sin(//,0/2)  —  0  but  from  (??),  we  then  find  n  —  0 
which  is  not  possible.  The  second  possibility  for  satisfying  Eq.  (??)  is  to 
have 

C2h  sin(//,0)  +  2 n  —  0.  (48) 

With  (??),  we  eliminate  Ch  in  (??)  and  obtain,  after  many  simplifications, 
and  for  /j,  /  0,  the  following  equation  for  p, 


sin(//.0/2)  = 
1  cos(p©/2) 


(49) 


which  we  recognize  as  Eq.  (??)  for  Bi  —  B-2.  We  conclude  that  p  = 
B-y  —  B2.  The  solution  of  the  linearized  problem  is  then  of  the  form  A  ~ 
(. Ri  +  cexp(ip^))  exp  [i  {Biu  +  dexp(ip;/))] which  we  may  rewrite  as 


A  ~  exp(fi?i^)  +  (c  +  w/i?i)  exp(fi?2^)  (50) 


for  small  c  and  d.  Thus  the  secondary  bifurcating  solution  is  close  to  a  linear 
combination  of  the  two  periodic  modes  at  the  double  Hopf  bifurcation  point. 


10  Figure  captions 

Figure  1.  Semiconductor  laser  subject  to  an  opto-electronic  feedback.  The 
figure  illustrates  the  opto-electronic  device  used  by  Saboureau  et  al  [?].  The 
feedback  operates  on  the  pump  of  the  laser  by  using  part  of  the  output  light 
which  is  injected  into  a  photodctcctor  connected  to  the  pump.  The  delay  of 
the  feedback  is  controlled  by  changing  the  length  of  the  optical  path. 

Figure  2.  Bifurcation  diagram  of  the  periodic  solutions.  Figure  2a  and 
Figure  2b  represent  the  amplitude  R  and  the  frequency  B  as  functions  of 
the  feedback  rate  C,  respectively.  Squares  denote  Hopf  bifurcation  points. 
The  values  of  the  parameters  arc  P  —  0.5,  T  —  1000  and  r  =  2000  (implying 
9-  63.25  and  0  =  4). 

Figure  3.  Hopf  bifurcation  points.  Figure  3a  and  Figure  3b  represent  the 
Hopf  bifurcation  point  Ch  and  the  Hopf  frequency  BH  as  functions  of  the 
delay  0,  respectively.  The  values  of  the  fixed  parameters  are  P  —  0.5  and 
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T  —  1000.  The  squares  denoted  the  Hopf  bifurcation  points  for  r  —  2000 
(0  —  63.25)  as  shown  in  Figure  2. 

Figure  4.  Primary  and  secondary  Hopf  bifurcation  lines.  The  Hopf  bifurca¬ 
tion  lines  C  —  Cp  (full  lines)  and  the  quasiperiodic  bifurcation  line  C  —  Cqj> 
(dashed  line)  are  represented  as  functions  of  the  scaled  delay  6.  The  points 
Pi  and  P4  are  double  Hopf  bifurcation  points.  The  point  P2  corresponds  to 
a  change  of  direction  of  the  Hopf  bifurcation  point  (supercritical  for  9  <  0p2 
and  subcritical  for  9  >  0p2).  At  the  point  P3,  a  Hopf  bifurcation  point  and  a 
quasiperiodic  bifurcation  point  may  have  the  same  value  of  the  bifurcation 
parameter  C  but  they  exhibit  different  amplitudes.  The  values  of  the  fixed 
parameters  are  P  —  0.5,  T  —  1000. 

Figure  5.  Qualitative  bifurcation  diagrams.  We  represents  the  amplitude  of 
the  rapid  oscillations  (R)  as  a  function  of  the  feedback  strength  (C).  The 
horizontal  line  represents  the  laser  steady  state.  It  changes  stability  at  a 
Hopf  bifurcation  point  which  leads  to  either  supercritical  oscillations  (Figure 
5a)  or  subcritical  oscillations  (Figure  5b  and  5c).  In  all  figures,  quasiperiodic 
oscillations  appear  as  the  result  of  a  secondary  bifurcation. 

Figure  6.  Numerical  bifurcation  diagrams.  Figure  6a  and  6b  represent  the 
bifurcation  diagrams  of  the  periodic  and  quasiperiodic  oscillations  obtained 
from  the  laser  equations  (??)  and  (??)  and  from  the  amplitude  equations 
(??)  and  (??),  respectively.  From  (??),  we  note  that  y  ~  2eAsin(.s  +  (f>)  and 
we  represent  max(y)  —  2 eR  in  Figure  6b. The  values  of  the  parameters  are  the 
same  as  in  Figure  2.  The  two  main  branches  of  periodic  solutions  and  their 
successive  bifurcations  have  been  obtained  by  changing  r/  by  successive  steps 
first  forward  and  then  backward.  The  dotted  lines  in  Figure  6b  correspond 
to  the  exact  expressions  of  the  branches  of  periodic  solutions  given  by  (??) 
and  (??). 
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